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ABSTRACT 

Imaging data from the Sloan Digital Sky Survey are used to characterize the population of galaxies in 
groups and clusters detected with the MaxBCG algorithm. We investigate the dependence of Brightest 
Cluster Galaxy (BCG) luminosity, and the distributions of satellite galaxy luminosity and satellite 
color, on cluster properties over the redshift range 0.1 < z < 0.3. The size of the dataset allows us to 
make measurements in many bins of cluster richness, radius and redshift. We find that, within r2oo of 
clusters with mass above 3 x 10^^h~^ Mq, the luminosity function of both red and blue satellites is only 
weakly dependent on richness. We further find that the shape of the satellite luminosity function does 
not depend on cluster-centric distance for magnitudes brighter than °-^^Mi - Slogip/i = —19. However, 
the mix of faint red and blue galaxies changes dramatically. The satellite red fraction is dependent 
on cluster-centric distance, galaxy luminosity and cluster mass, and also increases by ~5% between 
redshifts 0.28 and 0.2, independent of richness. We find that BCG luminosity is tightly correlated 
with cluster richness, scaling as Lbcg ~ -^2601 ^ Gaussian distribution at fixed richness, 

with aiogL ^ 0.17 for massive clusters. The ratios of BCG luminosity to total cluster luminosity 
and characteristic satellite luminosity scale strongly with cluster richness: in richer systems, BCGs 
contribute a smaller fraction of the total light, but are brighter compared to typical satellites. This 
study demonstrates the power of cross-correlation techniques for measuring galaxy populations in 
purely photometric data. 

Subject headings: galaxies: clusters — galaxies: evolution — galaxies: halos — cosmology: observa- 
tions 



1. INTRODUCTION 

Clusters of galaxies are important systems for study- 
ing both galaxy evolution and cosmology. Used as lab- 
oratories with well-defined environments, these massive 
objects are a tool for investigating processes that influ- 
ence galaxies' physical characteristics. Used as tracers 
of the underlying mass distribution, they are a tool for 
investigating the evolution of structure and the nature of 
dark energy in the Universe. These objectives are closely 
linked: cosmological studies require accurate knowledge 
of cluster selection and redshift- and mass-observable re- 
lations, and these facts are directly related the evolution- 
ary properties of the cluster galaxy population. 

In the context of galaxy evolution, the high-density en- 
vironment of galaxy clusters is a particularly interesting 
place to examine the galaxy population. Several studies 
have suggested substantial galaxy transformation in such 
environments. Galaxy morphology, star-formation rate, 
and luminosity have long been known to depend on clus- 
ter properties and to depart signiflcantly from the cos- 
mological average (e.ff., [Hubble 1926; Abcll 196^: lOemleJ 
1 1 9 741 : [Dressleij[l980l ) . Historically, the cluster galaxy con- 
tent has been quantified by the luminosity function (LF) 
and type fraction (such as the late- type or blue fraction) . 
Measurement of these quantities as a function of cluster 
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mass, redshift, and distance from the cluster center pro- 
vides insight into the underlying physical mechanisms 
responsible for these trends. 

While the LF is primarily a decreasing function of lu- 
minosity, galaxies are bimodal in color and spectral type, 
with red, early-type galaxies displaying little ongoing 
star formation, and blue, late-typegalaxics cxhibitin; 
signs of recen t star f ormation ('e.(/., [Strateva ot al. 200jj 
Baldrv et all |2004 iBell et all l2004r iMenanteau et' at] 
2006t iBlanton et al.1 l2005a[ ). This bimodality was in 



place by z ~ 1 at the latest and ma y provide a s i gnpos t 
of galaxy transformation {e.g., iFaber et al.l l2007f) . 
Although this bimodality persists in all environments, 
the fraction of galaxies in each class (a.k.a. the red 
and blue fractions) changes systematically with local 
density; this trend is the so-called morphology-density 
relationship ('Ocmler"l974l: lDresslei1ll980l : [Dressier et al.l 
[l997: Smith et al. 2005). In recent large galaxy surveys 
this work has been extended to show that a wide 
range of galaxy properties, including morphology, 
star-formation rate, and color, dep end on l ocal density 
(iGomez et al.l 120031: iBaloeh et al.l [2004allbl: iHogg et al.l 



20041: iKauffmann et alj 120041: iTovmassian et al 
Blanton et al.l l2005at IChristlein fc ZabludoS" 



Croton_ et all l2005t iRoias et all l2005l: ICooper et all 
2006; Ma ndelbaum et al.ll2006l : ICooper et al.ll2007l ). 

The advent of large galaxy surveys has made it possi- 
ble to place observational constraints on both the type 
fraction and the LF function of cluster mass (the 
conditional luminosity function, CLE), and to further in- 
vestigate how these quantities d epend on other va riables. 
In the Sloan Digita l Sky Survev (lYork et al.ll2000l . SDSS), 



iGoto et all (|2002t l and iHansen et al.l (|2005f l examined 



the LF as a function of cluster richness and cluster- 
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centric distance for systems found in the S PSS Early 
Data Release, and IWeinmann et al.1 ()2006bl ) measured 
the CLF measured from a group catalog derived from 
the spectroscopic sample of SPSS DR2. Us i ng the 2dF 
Galaxy Redshift Surve y, iDe Propris et al.l (|2003f ) and 
iRobotham et al.l (|2006( l compared the LF in high- and 
low-mass systems; a similar study was performed wit h 
a sample of 93 X-ray selected clusters ()Lin et al.ll2004[ ). 
The dependence of t he LF on galaxy co lor has been 
recently investigated (iPopesso c t aT '20 05D as has the 
galax y type fraction (|Goto et al.i,2003. : De Propris et al.l 
|2004| ) for clusters in these large surveys. 

The type fraction of galaxies in clusters depends 
on both cluster richness and redshift: the frac- 
tion of star-forming galaxies at fixed local density 
is larger at higher redshif t, an effect known as 
the Butcher-Oemler effect ([Butcher fc Oemleij Il978l 
[1981 . This effect is now well-documented, if 
not entirely well-explained, in clusters over a wide 



19911 iBTOTigh et al I [200l: IDe Grandi et all [2003 : 
Brough et all 120051 : iLoh fc Strausd I2006D . includ ing to 



or its 


converse, the red fraction (iRakos fc Schombert 


1995 




Margoniner & de Carvalhol 120001: lEUingson et al. 


2001 




Kodama fc Bower l2001l Margoniner et al. 2001: 


De Prooris et al.ll2004lMartmez et al.ll2006l: iGerke et al. 


20071). 


Other indicators of galaxy state, including 



galaxy morphology and e missio n line strength, also 
show a trend with r e dshift (lAUin gton-S mith et aDll993l; 
Oemleretall 119971: iBalogh et al.i .19971: ICouch etal.1 



1998L Ivan Dokkum et al 
Lubin et all l2002l: 



|20fld; iFasano et all I200C 

I: [Goto et al.l [20031: iTrcu ct al.' '200l 
Wilman et all 120051: iPoggianti et aLl [2006 : Dcsai ct ^ 
van der Wei et al.ll2007f ). However, few samples to 



date have had both large numbers of systems and well- 
understood mass proxies, so the dynamical range and 
mass resolution of these previous works has been some- 
what limited. 

Another characteristic of galaxy clusters is the pres- 
ence of a highly-luminous galaxy near the cluster center 
— the Brightest Cluster Galaxy (BCG). In addition 
to being extraordinarily luminous, BCGs differ in a 
number of ways from other clust er members: they tend 
to have exte n ded light p r ofiles ( Matthews et al.l 
Tonrvl Il987t ISchombertI Il988l : iGonzalez et al 



1964 



200C 



2003( 1 ■ larger size at fixed lu minosity than other early 



types ([Bernardi et al.l l2007l and references therein) 
and may contain a larger fraction of dark matter 
than typical galaxies (e.g ., iMandelbaum et all l2006l : 
Ivon der Linden et al.l 1 20071 ) . Also, while the t raditional 
fitting function to the LF of ISchechteil (|1976f ) provides 
a good fit for satellite galaxies, BCGs follow a different 
distri bution, causing the so-called "bright end bump" 
(e.g.. [Hansen et al]|2005f ). This difference between the 
BCG and other satellites in a cluster is also mani- 
fested in the luminosity gap statistic, the difference in 
luminosity between the BCG and the next brightest 
cluster member, which may be indicative of the special 
accretion history of B C Gs (lOstriker fc Tremaine Il975l : 
iTremaine fc Richstond Il977t [Loh fc StraussI l2006( ). 
BCGs are also distinct from other galaxies with sim- 
ilar mass that are not at the cente r of c luster-sized 
potential wells (|von der Linden et al.l l2007f ). Indeed, 
the properties of BCGs seem to be closely linked to 
properties of the ir hos t clusters (ISandage fc Hardvlll973l: 
ISchneider eral] 119831 : ISchombertI Il988i : lEdge fc Steward 



the masses of the i r par e nt halos ( [ Lin fc MohH |2004| : 
IMandelbaum et all l2006l : I Zheng et all l2007f) . The 
outer light profile of BCGs often merges smoothly 
with the diffuse intra-cluster light (ICL), suggesting 
again a coupling betw een the BCG and the host halo 
([Gonzalez et al.ll2005[ ). Generally t he BCG-HCL light 
is closely linked wi t h cluster mass (IZibetti et al.[ 120051: 
Conrov et all l2007l : iPurcell et all l2007l : IGonzalez et all 



2007ft . As BCGs have properties different from those of 



the rest of the cluster members, BCGs and satellites are 
often analyzed separately, and we follow this convention 
here. 

It is expected that the properties of cluster galax- 
ies are closely tied to the merging and accretion his- 
tory of their parent dark matter halos. Current mod- 
els based on Cold Dark Matter (CDM) suggest that 
the stars in BCGs were formed in dense peaks quite 
early but that the BCGs were assembled in a sc- 
ries of galaxy merging e vents that continue unti l rela- 



Dubinskil Il998l: iGao et al. 


2004: iBovlan-Kolchin et all 


20061: IDe Lucia fc BlaizotI 


20071). Satellite galaxies in 



clusters are now generally understood to be hosted by 
smaller dark matter halos that have merged into the par- 
ent halo. Several features of the observed cluster galaxy 
populations can be understood b ased on the assembly 
history of the parent halo {e.g., IPoggianti et al.l l2006l : 
iTasitsiomi et al.ll2007( ). One such characterization is the 
core distinction between central galaxies and satellites: 
the central concentrations of mass and light in massive 
halos continue to build up while the growth of satellite 
systems is halted upon accretion. BCG luminosity, and 
the correlations between this luminosity and both clus- 
ter mass and satellite properties, can thus provide insight 
into the assembly histories of clusters. 

The merger history of dark matter halos alone is not 
enough to account for the bimodality in galaxy properties 
and its dependence on redshift and environment. Several 
physical processes have been proposed to transform star- 
forming galaxies into the typical cluster galaxies on the 
red sequence. Although the relative strengths of these 
processes are still hotly debated, a consensus is emerg- 
ing that the main transformation mechanisms are related 
to the mass of the host halo and whether (and for how 
long) the galaxy has been a satellite within a larger sys- 
tem. Among the suggested transformation mechanisms, 
some are expected to be most effective in rich clus ters, 
such as ram-pressure stripping ([Gunn fc GottJfl97g)^ in- 
teraction with the cluster potential ([Byr d & Va ltonenI 
19901) and high-ve locity close encounters ("harassment;" 
Moore et al.lll996[ ). However, studies of very poor sys- 



tems have shown that the environmental dependence 
of galaxy properti es is not limited to the richest ob- 
jects ( e .q.,[Z abludo ff fc Mulchaevlll998l : IWeinmann et all 
[2006at [Gerice et all 120071 ). Processes that can operate 
efficiently in low velocity dispersion systems therefore 
must also play a role in shaping the galaxy po pula- 
tion: e.g., galaxy mergers ([Toomre fc Toomrel[l97l and 
"strangulation," a cutoff to gas acc retion onto galaxy 
disks by stripping or AGN feedback (iLarson et al.lll980l : 
IBalogh et aL[|2000l : ICroton et al.ll2006D . It is likely that 
some combination of these effects is at work. Distinguish- 
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ing their relative significance requires precisely quantify- 
ing cluster galaxy properties over a wide range of masses 
and as a function of cluster-centric distance. These data 
will provide information on both the assembly histories 
of clusters and on the physical mechanisms that trigger 
and quench star formation. 

Models of galaxy evolution in a cosmological context 
most readily predict galaxy properties as a function of 
halo mass rather than cluster observables. In order to 
make these comparisons, a reliable mass-obscrvablc re- 
lationship is a prerequisite. In addition, there is consen- 
sus that the luminosity function and type fraction both 
depend on a number of variables, complicating detailed 
comparison between clusters of different mass. For ex- 
ample, the cluster galaxy LF depends on cluster-centric 
distance, and the size of the bound regions of clusters 
scales with mass. In order to make physically meaningful 
comparisons between LFs of different mass clusters, an 
aperture scaled to the bound region is therefore prefer- 
able to a fixed metric aperture. With recent extensive 
surveys providing well-calibrated cluster catalogs span- 
ning a wide range in mass, it is possible to examine in 
detail the dependence of the cluster galaxy population 
on several cluster and galaxy properties simultaneously. 

Large, homogeneous photometric surveys such as the 
SDSS provide rich data with which to characterize the 
cluster galaxy population. These data have been used 
to define large, robust, clean samples of galaxy clus- 
ters with accurate photometric redshifts. These samples 
are sizable enough to split on several variables allow- 
ing detailed statistical exploration of the galaxy popula- 
tions in clusters. Currently, the largest sample of clus- 
ters available is the MaxBCG catalog fro m the Sloan 
Digital Sky Survey l|Koester et all l2007b| ). The selec- 
tion effect s of the cluste r-finding algorithm are well un- 
derstood (jKoester et al.. 2007a; Rozo et al. 2007a), and 
there are a number of studies ex ploring the mass- 
richness relat ionship for th ese objects (iRozo erani2007bl: 
Becker et aL|[200 7: Sheldon 'eralll2007al : iJohnston et ail 
2007bt iRvkoff et al.l [20071). In this work, we con- 
vert cluster richness to cluster mass using weak lens- 
ing measurements of the MaxBCG mass-richnes s rela- 
tion (jSheldon et al.ll2007at [Johnston et al.ll2007b[ ). The 
quality and quantity of these data allow for detailed mea- 
surements of a variety of cluster galaxy properties as a 
function of cluster mass and radius. For satellite galaxies 
we then measure the luminosity function of all, red, and 
blue satellites conditional on both mass and cluster ra- 
dius, and investigate the dependence of the red fraction 
of satellites on cluster mass, redshift, galaxy luminosity 
and distance from cluster center; for BCGs wc quantify 
the dependence on cluster mass of both the BCG lumi- 
nosity and the relationship between the BCG luminosity 
and satellite galaxy luminosities. Although this sample 
of clusters extends only to z = 0.3, these objects pro- 
vide a valuable low-redshift baseline with which higher 
redshift samples may be compared. 

In this paper we use a statistical background- 
subtraction technique to measure mean galaxy prop- 
erties over a wide range of color and luminosity in 
MaxBCG clusters. We average the signal from many 
clusters binned by cluster properties, and statistically 
subtract the contribution from random galaxies along 
the line of sight. This method of cross-correlating clus- 



ters with the galaxy population provides very precise 
statistical measurements, and allows us to study blue 
and low luminosity galaxies that are indistinguishable 
from the background in individual clusters. We test 
the background-correction algorithm by running the full 
analysis on realistic mock catalogs, and find that we 
are able to robustly recover 3D cluster properties us- 
ing these methods. The statistical techniques presented 
here for background correction of photometric data will 
be directly applicable to future large multi-band imag- 
ing programs, including the Dark Energy Survey (DES^ , 
lAljbott fc The PES Collaboration! I2OO5I) and t he Large 
Synoptic Survey Telescope (LSST6. [T7i^l2002D and the 
Panoramic Surv ey Telescope fc Ra pid Response System 
fPan-STARRS^ iKaiser et al.ll2002[ ) and will be essential 
for leveraging these data to provide the desired insights 
into cosmology and galaxy evolution. 

The paper is organized as follows: in §[2]we describe the 
SDSS and simulation data used; we present the stacking 
and background-correction method in § [S] Our primary 
results are given in § UJ § 14.11 presents the luminosity 
and color characteristics of the satellite population, while 
§ 14.21 discusses the BCG population. A summary and 
discussion of the implications of the results is given in § 

El 

The notation used for cluster-related variables in pre- 
vious MaxBCG work includes defining 7V200 and L200 as 
the counts and i-band luminosity of red-sequence galax- 
ies within the measurement aperture of the cluster finder 
and with L > 0.4L*. We note that as the aperture for 
cluster finding was determined with a previous definition 
of richness, it is not strictly the true value of r2oo for 
these systems (in fact, it is larger), and only red galax- 
ies arc included in these definitions. In this work we 
will refer to the total excess luminosity associated with 
the light from galaxies of all colors above a luminosity 
threshold and within the measured r2oo of these systems 
as ^200- In addition, wc follow the standard convention 
of using R to to denote projected, 2D radii and r to refer 
to deprojected, 3D radii. Where necessary for comput- 
ing distances, we assume a flat, LCDM cosmology with 
Ho = 100/1 km s^^ Mpc^^, h = 0.7, and matter density 

^rn - 0.3. 

2. DATA 

We use the photometric data of the Sloan Digital 
Sky Survey (SDSS) Data Release 4. The SDSS data 
were o btained using a specially-designed 2.5-m tele- 
scope (jGunn et al.l l2006( ) , operated in drift-scan mode 
in five bandpasses (u, q,r,i, z) to a lirn i ting magni- 
tude of r < 22.5 ^Fukugita et all Il996l: iGunn et all 
19981 iLupton et all lin preparationt iHogg et al.l l200lF 
Smith et al.ll2002l ). The survey covers much of the North 
Galactic Cap and a small, repeatedly-scanned region in 
the South. Apparent magnitudes are corrected fo r Galac- 
tic extinction using the maps of iSchlegel et"an (fl99l . 
Photometric errors at bright m agnitudes are < 3 %, lim- 
ited by systematic uncertainty (jivezic et al.ll200^ ; astro- 
metric er rors are typically smaller than 50 mas per co- 
ordinate (jPier et al.ll2003f) . Further details of the SDSS 

^ http : //www ■ darkenergysurvey ■ org | 
^ http: //www. Isst . orgi 
http : //pan-Starrs . if a . hawall ■ edu | 
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d ata may be found in the most rece nt data release paper 
of lAdelman-McCarthv et all (|2007| ). 

With these data, we define a catalog of galaxies, then 
run our cluster-finding algorithm on this list to produce a 
catalog of galaxy clusters. The set of clusters and galax- 
ies are then jointly analyzed to determine the binned 
correlation function measurements used in this work to 
describe the galaxy population associated with clusters. 
In this section, we describe the galaxy and cluster cata- 
logs. 

2.1. Galaxy Sample 

Our photometric galaxy catalog was generated by ap- 
plying the Bayesian star- g alaxy separation method de- 
veloped in lScranton et all (l2002f) to the full set of SDSS 
data. The primary source of confusion in star-galaxy sep- 
aration at faint magnitudes is shot noise, which causes 
stars to scatter out of the stellar locus and galaxies to 
scatter into the stellar locus. The IScranton et al.l ()2002f ) 
method uses knowledge about the underlying size distri- 
bution of objects as a function of apparent magnitude 
and seeing to assign a probability of being a galaxy to 
each object. We examined objects from the repeatedly 
scanned SDSS southern stripe and co-added the fluxes 
from different exposures at the catalog level to provide 
more precise sizes and magnitudes. Using regions with 
at least 20 scans, we selected objects on single exposures 
with typical observing conditions, but characterized their 
size distribution as a function of magnitude and seeing 
with the co-add. From these distributions we generated 
a galaxy probability for every object in the survey. The 
resulting distribution is highly peaked at probabilities 
and 1, such that prob > 0.8 results in a highly pure sam- 
ple of galaxies. We then searched the resultant catalog 
of galaxies for systems that are likely to be clusters of 
galaxies using the Ma xBCG algorithm desc ribed briefly 
below and in detail in iKoester et al.l (|2007b[ ) . 

For characterizing the galaxy population in clusters, 
we restrict the galaxy sample to be volume and magni- 
tude limited for z < 0.3 by including only those galax- 
ies brighter than an appropriate threshold. Expressed 
as a luminosity K-corrected to redshift z = 0.25 (see § 
13. 4p . this lower i-band luminosity limit is lO^'^/i~^L0. 
At z = 0.3, this value corresponds to an apparent mag- 
nitude limit oi i < 21.3, and color limits of '^'^^{g — r) 
< 2 and 0-25 - i) < l. AH magnitudes are SDSS model 
magnitudes. Uncertainties on the measurements are suf- 
ficiently small so that the typical r — i signal-to-noise is 
greater than 30 until beyond redshift 0.3. Galaxy col- 
ors are the refore robus tly measured (ISmail et al.l 0,9951 : 
iBernstein e t al. 2002; Benitez et al.l |2004[) overthe full 
redshift range considered here. 

2.2. Cluster Sample 

Cluste rs were ident i fied us ing the MaxBCG cluster 
finder of lKoester et al.l ()2007bD . This algorithm identifies 
clusters by the presence of a BCG and a red sequence, 
and provides an excellent photometric redshift for each 
system. The scatter in redshift is Az < 0.02 for our 
full richness and redshift range and Az ^ 0.015 for A^200 
> 10. The cluster center is defined to be the position 
of the BCG. The cluster richness, iV2oo, is defined as the 
number of red galaxies with rest frame i-band luminosity 
L > 0.4L, in an aperture that scales with r2oo, and is an 



excellent mass proxy, as discussed below. The catalog of 
system s with iV2oo > 10 was presented in IKoester et al.l 
()2007al) : in this work we use systems with A^200 > 3. 

The MaxBCG cluster finder was trained using known 
BCGs and luminous red galaxies, and as we will con- 
sider whether the resulting BCG population is being con- 
strained by these priors, we review the procedure here. 
Specifically, when constructing the likelihood of a galaxy 
to be a BCG, the algorithm places priors on observed g—r 
and r — i colors and i-band magnitude. The model BCG 
color-redshift and magnitude-redshift relationships are 
taken from the distribution of bright LRG galaxies, then 
refined based on examination of 100 visually identified 
BCGs in rich clusters. The BCG likehhood, £bcg, is 
specified to be 

£bcg{z) = Gff,«(z)Gfi;;«(z)e-«™— )/-)^ (1) 

(jKoester et al.ll20"07bL Equation 11), where the width CTc 
of th e prior on BCG i-band magnitude is 0.3mag, taken 
from iLoh fc StraussI (|2006( ). In addition, BCGs are re- 
quired to be brighter than 0.4L, (6.9 x lO^/i'^L© at 
z = 0.25). The distributions G^5/^ and Gf^^^ are nar- 
row Gaussians in color; the width of each includes both 
the intrinsic width of the E/SO ridgeline (0.05 for g — r 
and 0.06 for r — i) and the error in the measured color 
for the galaxy in question. In practice, the narrow color 
priors provide a much stronger constraint in identifying 
likely BCGs than does the magnitude prior, and for the 
majority of cases, the BCG is clearly defined by the color 
criteria. The BCG magnitude likelihood has an appre- 
ciable effect only for the cases where the choice of BCG 
is ambiguous. 

The selection function of the MaxBC G alg orithm was 
investig ated by IKoester et al.l (|2007b( l and iRozo et all 
([2007^, who showed that this sample of clusters is over 
90% complete and pure for N200 > 10. For N200 > 20 
sytems, the completeness and purity rates are 97% or 
better. The completeness and purity were quantified 
by comp aring clusters identifie d in realistic mock galaxy 
catalogs (jWechsler et al.ir2007( ) to their host halos. We 
did not wish to use unresolved halos in our comparisons, 
so for lower richness values (A''2oo < 10), the complete- 
ness and purity are less well quantified. It is likely that 
the lowest-richness systems are less pure and complete, 
not only because of statistical considerations but also 
because the lower richness systems, with just a few red 
galaxies in close proximity, are not necessarily represen- 
tative of all systems in the equivalent mass range. Higher 
resolution simulations will be necessary to quantify this 
selection in detail; however, we can learn about the statis- 
tical properties of these systems even if they are a biased 
population compared to dark matter halos. We thus in- 
clude them in this study, but for the purpose of fitting 
relationships of galaxy properties as a function of cluster 
richness or mass, we restrict the sample to iV200 > 10, 
where the catalog is known to be both highly complete 
and pure. 

We use clusters in the redshift range 0.1 < z < 0.3; 
there are 165,597 objects identified with A^200 > 3, of 
which 13,823 have iV200 > 10. 

2.3. Average Galaxy Properties at z = 0.25 

Below, we compare our derived cluster galaxy pop- 
ulation statistics to similar statistics for the mean of 
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0.085 Mpc 



r = 0.034 Mpc 



r = 0.049 Mpc 



0.070 Mpc 



r = 0.10 Mpc 




10.0 10.5 11.0 11.5 10.0 10.5 11.0 11.5 10.0 10.5 11.0 11.5 

log(Lo.25i/L0h^) 

Fig. 1. — Joint distribution of '^'^^(g — r) and '' ■^^i-band luminosity density as a function of projected separation from cluster centers for 
an example richness bin, 12 < Af200 < 17. Each panel corresponds to a different radial bin as indicated. The one-dimensional distributions 
for luminosity and color arc shown along the x- and j/-axes respectively. The luminosity distribution is expressed as log of the number 
density as a function of log luminosity; the color distribution is given as the linear density as a function of color. At small separations, 
red galaxies dominate, while on large scales there is a bivariate color distribution similar to the global distribution. A smaller fraction of 
galaxies are highly luminous near cluster centers. 



galaxies in all environments, which we will refer to as 
the cosmological, or universal, average. To calculate 
these average galaxy statistics, we use the New York 
University Value-Ad ded Galaxy Catalog (NYU-VAGG 
iBlanton et al.|[2005bD of the SDSS spectroscopic sample, 
whose median redshift is z ~ 0.1, evolved to z = 0.25. 
The details of constructing this sample are discussed in 
iSheldon et all (|2007bf ). By applying the same luminos- 
ity and color cuts to this sample as we do to any given 
set of cluster galaxies {e.g., to determine the density of 
red, bright galaxies), wc estimate the universal average 
values of a comparable "field" sample. We do not un- 
dertake a detailed study of this population here, but just 
present the average values for general comparison with 
the recovered cluster values. For reference, from fitting 
a Schechter function to this sample, wc find the global 
value of evolved to z = 0.25 to be 1.7 x 10^°h~'^LQ. 

3. METHOD 

In this section, wc describe how wc characterize the 
radial, luminosity, and color distribution of galaxies that 
are associated with clusters. In 13.11 we review the esti- 
mator used to determine the distribution of cluster galax- 
ies; §§ 13.21 through 13.51 discuss the technical details of 
characterizing the survey geometry and redshift distribu- 
tions, performing K-corrections, and determining M200 
and r-zoQ . Wc describe a check to the background correc- 
tion with a mock galaxy catalog in § 13.61 



We can only reliably determine cluster membership for 
individual clusters in this purely photometric sample for 
red, relatively high luminosity galaxies. We nonetheless 
can accurately characterize the full cluster galaxy popu- 
lation by statistically correcting for galaxies aligned by 
chance along the line of sight. We perform this correc- 
tion not on individual clusters, but for ensembles binned 
together based on their observable properties, such as 
richness or redshift. We do not recover the properties of 
individual clusters, but rather properties of the average 
cluster. In this sense we measure the cross-correlation 
between clusters and the galaxy population. 

If the average stacked cluster is spherically symmetric, 
we can invert the projected densities to recover the 3D 
profile. This assumption is true for a homogeneous and 
isotropoic universe, as long as the cluster finder does not 
introduce anisotropics: for example, choosing preferen- 
tially systems elongated along the line of sight. 

To recover the K-corrected luminosities of the average 
galaxy population, we K-correct galaxies along both the 
cluster line-of-sight and the random lines-of-sight to the 
redshift of the cluster, regardless of their true redshifts. 
After background subtraction, galaxies at redshifts dif- 
ferent from the cluster are removed from the counts, and 
so this inaccuracy does not persist in our final measure- 
ments. 

A correlation function-based approach is straightfor- 
ward to implement and interpret, and takes advantage of 
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the full dataset for background correction. This method 
is in contrast to using a "local" background correction 
estimate from, for example, counting galaxies in a large 
annulus around each cluster, as has been done in studies 
of more limited area. As there is mass and light corre- 
lated with clusters out to many times the virial radius, 
a local background estimate is likely to result in a bi- 
ased estimate of the true galaxy population in clusters. 
iMasiedi et al.l (|2006[ ) implemented the use of a correla- 
tion function-based estimator to investigate galaxies cor- 
related with the SDSS LRG sample, and here we extend 
this technique to examine galaxies associated MaxBCG 
clusters. This estimator is a powerful tool, but requires a 
significant amount of computation for the large number 
(here 2.7 billion) of pairs that must be considered. 

3.1. Cluster- galaxy Correlations 

We follow the method of lMasiedi et al.l (|2006D to esti- 
mate the mean number density and luminosity density 
of galaxies associated with clusters. This method es- 
sentially calculates a correlation function with units of 
density; it includes corrections for random pairs along 
the line of sight as well as pairs missed due to edges and 
holes. In the description of the estimator below, we dis- 
cuss the necessary terms for calculating the luminosity 
density. The number density measurement is simply a 
special case where we weight each galaxy by unity rather 
than by its luminosity. 

We define two samples: the primary sample, denoted 
p, and the secondary sample, denoted s, that are either 
in the real data or random locations Rx- For exam- 
ple, the counts of real data secondaries around real data 
primaries is denoted DpDs, while the counts of real data 
secondaries around random primaries is RpDs- For this 
study, the real data primaries are galaxy clusters with 
redshift estimates and the secondaries are the imaging 
sample of galaxies with no redshift information. The 
generation of random data is discussed in § 13.21 

The mean excess luminosity density of secondaries as- 
sociated with primaries is 



£ w 



DpDs 
DnRs 



RpDs 
RpRs 



(2) 



where £ is the mean luminosity density of the secondary 
sample, averaged over the redshift distribution of the 
primaries, and w is the projected correlation function . 
This quantity is the estimator from lMasiedi et ahl ()2006l ) 
where the weight of each primary-secondary pair is the 
luminosity of_the secondary. We have written the mea- 
surement as £w to illustrate that the estimator gives the 
mean density of the secondaries times the projected cor- 
relation function w. Only the excess density with respect 
to the mean can be measured with this technique. Using 
a weight of unity for each galaxy, rather than luminoisty, 
gives the mean number density. 

The first term in equation [2] estimates the total lumi- 
nosity density around clusters, including everything from 
the secondary imaging sample that is projected along the 
line of sight. The second term estimates the contribu- 
tion from random objects along the line of sight. Note, 
the same secondary may be counted around multiple pri- 
maries (or random primaries). 

The numerator of the first term, DpDg, is calculated 



DpD, 



J2p,s 



= {Lpa„.) + {Lt^r) = fA I {W + 1), 



(3) 

where the sum is over all pairs of primaries and secon- 
daries, weighted by the luminosity of the secondary. The 
secondary luminosity is calculated by K-corrccting each 
secondary galaxy's flux assuming it is at the same red- 
shift as the primary (see § 13.41 for details). The total 
luminosity is the sum over correlated pairs (Lpair) as 
well as random pairs along the line of sight [L^air)- -^y 
the definition of w, this number is the total luminosity 
per primary times w + 1. This term can be rewritten in 
terms of the luminosity density of the secondaries £ times 
the area probed A. Some fraction of the area searched 
around the primaries is empty of secondary galaxies due 
to survey edges and holes. The geometry factor / rep- 
resents the mean fraction of area around each primary 
actually covered by the secondary catalog. The geome- 
try factor is a function of pair separation, with a mean 
value close to unity on small scales but then dropping 
rapidly at large scales. Measurement of the survey ge- 
ometry is discussed in § 13.31 

The denominator of the first term calculates the factor 
f A, the actual area probed around the primaries. This 
term in the denominator corrects for the effects of edges 
and holes. Also, because the denominator has units of 
area, we recover a volume density rather than just the 
correlation fimction. This term is calculated as: 



DpRs 



pair 



E( dn\ dN_ 
p KdAJp dn 



fA. 



(4) 



The numerator is the pair counts between primaries and 
random secondaries, and the denominator is the expected 
density of pairs averaged over the redshift distribution of 
the primaries, times the number of primaries. The ratio 
is the actual mean area used around each primary fA. 

The second term in equation [2] accounts for the ran- 
dom pairs along the line of sight. The numerator and 
denominator of this term are calculated in the same way 
as the first term in cquation[2l but with randomly chosen 
locations as primaries distributed over the survey geom- 
etry. The redshifts are chosen such that the distribution 
of redshifts smoothed in bins of Az = 0.01 match that of 
the clusters. The numerator and denominator are 



RpDs 



y 



Ls 



R aR 
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and 
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dN 
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rA 



(5) 



(6) 



respectively. The ratio of these two terms, RpDg/ RpRs, 
calculates the mean density of the secondaries around 
random primaries after correcting for the survey geome- 
try- 

The density measured with this technique may be tab- 
ulated in various ways, typically as a function of pro- 
jected radial separation R. We tabulate in a cube that 
represents bins of separation R, luminosity ^ "^^Li, and 
color "■■^^(5 — r). This choice facilitates the study of 
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Fig. 2. — Mean fractional area searched relative to the area in 
the radial bin for one bin in cluster richness. Edges and holes dilute 
the true galaxy counts, biasing the measured density. This effect 
is negligible on small scales, but becomes important for large dis- 
tances from clusters, where a higher fraction of clusters are closer 
to the edge than the search radius. Due to the small area enclosed 
at small separations, the value is not well determined, but must ap- 
proach unity smoothly. We model this behavior with a polynomial 
constrained to unity on small scales. For larger scales no model is 
needed. 



the radial dependence of the luminosity function and 
the color-density relation. We make measurements for 
objects in the ranges 0.02 < R < 11.5h~^ Mpc (in 
18 bins), < °-^^{g - r) < 2 (in 20 bins) and 9.5 < 
logio(°^ i^/i©) < 11-7 (in 20 bins), where g - r and 
Li were K-corrected to the median redshift of the clus- 
ter sample, z = 0.25. Note that the smallest separation 
that we use is 20h~^ kpc, interior to which the light is 
dominated by the BCG. We retain information about the 
BCG of each cluster separately. An example of the result- 
ing distributions is shown in Figure [1] each panel shows 
the joint color-luminosity distribution in the richness bin 
12 < iV2oo < 17 for a different radial range, the mean of 
which is indicated in the legend of each panel. We save 
the cubes of luminosity density and number density data 
for every cluster, then bin the clusters into samples based 
on cluster observables. Wc evaluate the uncertainties for 
our results using jackknifc estimation from dividing the 
survey area into 12 roughly contiguous, roughly equal 
area sections. 

For each bin in richness, the projected radial profiles 
as a function of luminosity and color are inverted us- 
ing a standard Abel type integral to recover the three- 
dimensional profiles. The inversion relates the derivative 
of N{R), the projected number density profile, to n{r), 
the 3D density profile, as 



An(r) — n{r) — n 



dR 



-N'{R) 



(7) 



where An(r) is the excess density over the mean ran- 
dom background n. Inherent in this inversion is the as- 
sumption that N{R) is the linc-of-sight projection of a 
spherically symmetric n{r). In an isotropic universe this 



assumption is valid as long as cluster selection does not 
preferentially find systems oriented in a particular direc- 
tion with respect to the observer. 

We a pply the nonparametric method of lJohnston et al] 
(|2007af ) to deproject the number density profiles. Be- 
cause the maximum separation we use is lOh^^ Mpc, 
not infinity, we cannot constrain the 3D density at the 
outermost point, and the second to last point, at ~ 8h~^ 
Mpc, must be corrected. With a power law extrapolation 
of the profile, the correction is ~ 5% upward for that 
point. Deprojection is necessary to accurately recover 
the the red fraction and luminosity functions. Failure to 
properly deproject will result in artificially high faint-end 
LF sl o pes and artificially low red fractions (jValotto et al.l 
120011: iBarkhouse et al.l |2007[) . After deprojccting, we 
translate the bins of physical radial separation into bins 
of r/r2Qo using the measured r2oo-A^20o relationship (see 
§ 13. 5p . To make measurements within the same bins of 
r/r2oo for clusters of different richness, we use power-law 
interpolation between the fixed R values tabulated in the 
data cube. 

We construct the LF of clusters in a given richness 
range in the following manner. For each richness value, 
for a given bin in luminosity, we calculate the excess num- 
ber density of galaxies in the desired color and r/r2oo 
range. For the LF of a broader richness bin, we average 
the LFs of all clusters in that richness range. To measure 
the red fraction for clusters in a given richness range, we 
follow an analogous procedure: we find the number of 
red and all galaxies in the desired luminosity and r/r2oo 
range for each cluster richness, then take the average over 
all clusters within the broader richness bin. The number 
of red galaxies divided by the total number of galaxies 
(within the luminosity and r/r2oo limits) is the red frac- 
tion, /r. To measure the total luminosity of galaxies in 
a given r/r2oo, color, luminosity and richness range, we 
again follow the same procedure, but using the tabulated 
values of luminosity density, rather than number density. 

These measurements allow investigation of the distri- 
bution of excess-over-random galaxies around the BCGs 
used as MaxBCG centers. If a significant fraction of 
BCGs are miscentered relative to the density peak of the 
dark matter halos (whether because the true BCG does 
not lie at the deepest part of the potential well or because 
the wrong galaxy was identified as the BCG during clus- 
ter detec tion), the resulting we ak lensing profiles can be 
affected () Johnston et al.ll2007af l. as should the hght pro- 
files. In this work we make no miscentering correction, 
but remind the reader that we measure the distribution 
of galaxies around BCGs, not halo mass peaks. 

3.2. Random Catalogs 

To correct for the contribution of galaxies along the 
lines of sight to clusters, we generate a set of 15 million 
random points. These points are used in the RpD^, DpRg 
and RpRs terms from the estimator described in equa- 
tion [2l The random positions are distributed uniformly 
over the survey area using the window function described 
in § 13.31 The redshifts used for random points must sta- 
tistically match that of each cluster sample in order for 
the background subtraction to be accurate. To achieve 
the proper distribution, the random primaries were gen- 
erated with constant comoving density. For a given sub- 
sample of clusters wc re- weight the random primaries so 
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that the weighted redshift distribution matches that of 
each cluster subsample when binned with Az = 0.01. 

3.3. Survey Geometry 

We characterize the survey geometry using the 
SDSSPix code^. This code represents the survey using 
nearly equal area pixels, including edges and holes from 
missing fields and "bad" areas near bright stars. We re- 
move areas with ext inction greater than 0.2 magnitudes 
in the dust maps of ISchlegel et all {1998). This window 
function was used in the cluster finding and in defining 
the galaxy catalog for the cross-correlations. By includ- 
ing objects only from within the window, and generat- 
ing random catalogs in the same regions, we control and 
correct for edges and holes in the observed counts as de- 
scribed in § 13.11 The resulting area examined is 7398.23 
deg^. 

The denominator terms DpRg and RpRs in equation 
[2] correct for the survey edges and holes by measuring 
the actual area searched. An example DpRs is shown 
in Figure ^ generated for one richness bin. This quan- 
tity is expressed as the mean fractional area searched 
relative to the area in the bin. For small separations, 
edges and holes make little difference, so the fractional 
area searched is close to unity, but on larger scales edges 
become important. 

On very small scales the area probed in each bin is 
relatively small and the correction factor is not well con- 
strained. However, we know that the area correction 
factor must be close to unity, a fact that is clear from vi- 
sual inspection. To smooth the result, we fit a fifth order 
polynomial, constrained to be unity on small scales, to 
the fractional area as a function of the logarithmic sepa- 
ration. Due to the weighting, this approach results in a 
curve that approaches unity smoothly on small scales, yet 
matches intermediate separation points exactly. Points 
on larger scales are well-constrained and do not need 
smoothing. 

3.4. K-corrections 

We calculate d K-corrections using the template code 
kcorrect from iBlanton et all (j2003a[ ). This code is ac- 
curate but too slow to calculate the K-corrcctions for the 
over 10^ pairs found in the cross-correlations. To make 
the computation tractable, we computed K-corrcctions 
on a grid of colors in advance using galaxies from the 
SDSS Main sample as representative of all galaxy types. 
We computed the K-corrections for these galaxies on 
a grid of rcdshifts between and 0.3, the largest red- 
shift considered for clusters in this study. The mean K- 
correction on a 21x21x80 grid of observed 5— r, r— i, and z 
was saved. We interpolated in this cube when calculating 
the K-correction for a given galaxy. This interpolation 
makes the calculation computationally feasible for this 
study, but is still the bottleneck. 

To minimize uncertainties in the K-correction, we K- 
correct to the median redshift of the cluster sample, 
z = 0.25. All reported magnitudes are a djusted to this 
redshift , and arc noted as e.g.^ °-^'^Mi. iBlanton et al.l 
(|2003bf ) give a complete discussion of this bandpass shift- 
ing procedure. 

® |http: //lahmu.phyast .pltt . edu/'^scrantoii/SDSSPix/| 
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Fig. 3. — Cluster radius r200 as a function of cluster richness, 
measured with three different methods. The result from this paper, 
based on the distribution of galaxies, is shown as the points with 
error bars; the solid line is the best fit power law to systems with 
^200 > 10 and the dashed line extends this relationship to lower 
richness. These data are comp ared with results fr om velocity dis- 
persion profiles (hashed region: [Becker et al.||2007l1 and from weak 
lensing (shaded region; IJohnston et al.ll2007bir " The error bars on 
the galaxy data points are determined from jackknife resampling 
of the data. The width of velocity dispersion-based estimate cor- 
responds to an 18% uncertainty in mass; the width of the lensing 
region results from the reported 13% uncertainty in mass. 



3.5. r2oo and M; 



200 



There has been much recent work to quantify the 
mass-observable relation for stacked samples of the 
maxBCG cluster s, using mass est imates from clus- 
ter a bundance (iRozo et al.l l2007bf ). ve locity disper- 
sion (Becker et all I2007D. we ak lensing (jSheldon et al.l 



[2007i 



a; 



.Johnston et al. f l2007bD . ray measurements 

( Rvkoff et al.ll2007D . These various methods all result in 
a consistent mass-richness scaling; a detailed compari- 
son of the different mass estimators is in progress (Rozo 
et al. in preparation). These measurements also yield 
estimates of cluster size. 

To compare equivalent regions of clusters of various 
masses and thus various sizes, we scale all physical radii 
to r200j where r2oo is the threshold radius interior to 
which the mean mass density of a cluster is 200 times the 
critical mass dens ity of the Unive r se. The weak lensing- 
measured r2oo of I Johnston et all ()2007bf ) yields a size- 
richness relationship of 

(8) 



r2oo = 0.182;i-iMpc7V2°of 



which we use here. As we have measured 3D profiles 
around all clusters out to 8h~^ Mpc, we can construct 
3D radial profiles out to 5 x r2oo for even the most massive 
systems in our catalog. 

For compa rison, we also derive r2oo using the approach 
presented in [Hansen et al.l (|2005| ) for estimating r2oo as 
a function of A^200 using the space density of galaxies 
in clusters. In that work, the radial space density pro- 
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Fig. 4. — A test of the background-correction algorithm with the 
ADDGALS mock catalog. Each panel shows the difference between 
the underlying 3D distribution of luminosity and the background- 
correction algorithm's recovered distribution, in units of the un- 
certainty (T. The left column shows galaxies within 0.3xr2oo of 
cluster centers, while the right column shows galaxies in the ra- 
dial range 0.3 < r/r2oo !^ 0.5. The top panels use galaxies around 
clusters in a low richness bin (25 < N200 < 40), and the bottom 
panels use clusters in a high richness bin (60 < N200 < 100). In 
each panel the difference, in units of cr, is shown for the red fraction 
(purple asterisks), the density of all galaxies (black, filled circles) 
and of red galaxies (red, open diamonds). For magnitudes brighter 
than ^'^^Mi - 51ogio/i ~ -23, no galaxies are found, and so the red 
fraction is not calculated. In each case the background-correction 
algorithm accurately recovers the underlying distribution. 



file of cluster-correlated galaxies brighter than a luminos- 
ity threshold was compared with the global mean space 
density as d etermined from the g lobal SDSS luminos- 
ity function (|Blanton et al.|[2"003b[ ). evolved to z = 0.25, 
integrated to the same luminosity threshold. The ra- 
dius interior to which the mean space density of cluster 
galaxies reached 200/nM times the global mean density 
of galaxies was taken as r2oo ■ If galaxies were completely 
unbiased with respect to dark matter on all scales, r2oo 
measured with galaxies should exactly match r2oo mea- 
sured from the true mass distribution. The characteristic 
cluster-associated galaxy used has a luminosity of slightly 
sub-X*. On average in the Universe, a galaxy with this 
luminosity has a bias of close to unity. As our sample is 
specifically chosen to be galaxies associated with clusters, 
however, and clusters do not present the same environ- 
ment as the cosmological average, we cannot necessarily 
expect that the typical cluster galaxy thus also has bias 
close to one. Nonetheless, we show below that using these 
galaxies as a tracer of the total matter distribution docs 
result in a re asonable approximatio n of ^200, a finding 
also noted bv lGonzalez et al.] (|2007f ). 

Using the galaxy distribution-based method with the 
current data set, and comparing to the z = 0.25 field 
sample, we find for iV2oo > 10 that r2oo = (0.224 ± 
0.004)/i-i Mpc N^^J^°-°^. Figure El shows r2oo as es- 
timated from the galaxy distribution (points with er- 
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Fig. 5. — Luminosity functions of galaxies within r2oo for four 
example bins in cluster richness, showing the contribution from 
BCGs (open diamonds) and satellite galaxies (filled circles) to the 
total LF (solid lino). The figure illustrates one way in which the 
BCG population is different from the satellite distribution. Error 
bars are from jackknifc resampling, and are omitted on the total 
distribution (BCG-f satellites) for clarity. 



ror bars), the relationship from the velocity dispersion- 
based mass estimate (hashed region), and the lensing es- 
timate (shaded region). The width of the hashed region 
corresponds to a 18% uncert ainty in the mean mass- 
observable relationship from iBecker et al.l (j2007t ): the 
width of the shaded region corresponds to a 13% un- 
certainty in ma s s as re ported for the lensing analysis in 
iJohnston et^ (|2007bf ). The best fitting power law for 



-^200 ~ 10 for the galaxy-based measurement is plotted 
as the solid line, with a dashed line extending the trend 
to lower richness. The r2oo derived from the galaxy dis- 
tribution is remarkably consistent with these other ob- 
servations. The slight discrepancy in scaling gives hints 
about the way that galaxies populate halos of dark mat- 
ter, but further investigation of this interesting problem 
is beyond the intended scope of this work. 

The measurement of r2oo also enables an estimate of 
cluster mass, M2oo- For the present work, we present 
our results as a function of the direct observable, 
but use the r2oo^.^20o relationship from the weak lensing 
analysis to translate this observable into a mass estimate. 
The mass is calculated from r2oo via 
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200 



200pe(2)g7rr: 
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(9) 



where Pc{z) = 3i?^(z)/(87rG) is the critical density at 
epoch z, and the Hubble parameter is given by H^{z) = 
i?^[0^(l + z)3 + {1- Qra)] for a flat LCDM Universe. 
Fo r the choices of z = 0.25 and fi™ = 0.27 used 
in iJohnston et al.l jj2007bD , the M2oo^''20o conversion is 



Af2oo = 2.923 X 10 
fore adopt 



*/i-iMQ(r2oo//i"^Mpc)3. We there- 
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Fig. 6. — Color-luminosity distributions for satellite galaxies in clusters, in several example richness and radial ranges. Shading shows 
the relative density in color— luminosity space for each bin. Left: galaxies with 0.05 < r/r2oo ^ 0.3; right: galaxies in the range 3 < r/r2oo 
< 5. The dashed line indicates the (fixed) split used to separate galaxies into red and blue samples, given by equation 1111 

correction algorithm run on the mock. The mass Umit 
for resolved halos in the simulation, A/200 = 5 x 10^'^ 
h~'^MQ, allows a close to complete sample for a richness 
of A''2oo > 10, so we focus on results above that range. 

We find that the background-correction algorithm does 
well at reproducing the underlying distribution. For ex- 
ample, Figure[4]dcmonstrates that the difference between 
the underlying 3D values and the background-corrected 
values for several statistics are typically within ~ la. 
Shown is the comparison between the luminosity func- 
tions of all galaxies (filled circles) and red galaxies (di- 
amonds). Also shown is the comparison of red fraction 
values (asterisks). In all cases the difference is expressed 
in units of the uncertainty a. The left column is for 
galaxies within 0.3xr2oo, while the right column is for 
galaxies in the range 0.3 < r/r2oo < 0.5; the top panels 
are a low richness bin (25 < N200 < 40) and the bottom 
panels are a high richness bin (60 < A^20o < 100). For 
magnitudes brighter than ^-"^^Mi - 51ogio/i ^ -23, there 
are no galaxies found, so the red fraction is not calcu- 
lated. 



the best-fit value from I Johnston et all (|2007bf ). We note 
that these values are measured in the lensing data includ- 
ing the effects of miscentering, and so should reflect the 
underlying halo mass, rather than the mass correlated 
with cluster centers. 

3.6. Testing with Simulation Data 

To check that our methods of background correction 
and deprojection are reliable, we compare the SDSS data 
to results from a mock catalog from an N-body sim- 
ulation popu lated with galaxies u sing the ADDGALS 
method of of IWechsler et al.l ()2007f ). This catalog is de- 
rived from_the_N^od;Y_Hubble Volume light-cone simula- 
tion (jEvrard et al.l[2003 ). Particles in the simulation are 
assigned luminosities to match the luminosity-dependent 
two- point corre lation function measured in the SDSS 
bv Zchav et al. (20 05) an d th e global SDSS luminosity 



function ( Blanton et al.ll2003b[ ). Colors are assigned de- 
pending on local environment as defined by the distance 
to the 5th nearest neighbor, matching the photometric 
properties of SDSS galaxies with similar luminosities and 
local densities. The luminosity limit of the mock cata- 
log is L > 0.4L*, matched to the luminosity threshold 
used for identifying MaxBCG clusters. This mock cata- 
log successfully reproduces the width, location and evo- 
lution of the red ridgeline, making it ideal for testing and 
understanding the selection function of the MaxBCG al- 
gorithm. A detailed characterization of the MaxBCG 
selection function based on thi s cat alog has been under - 
taken bv lKoester et al.l (|2007bf ) and lRozo et"aI1 (|2007aD . 

We use this simulation to test the algorithms pre- 
sented for quantifying the galaxy population. We run the 
MaxBCG cluster finder on the mock catalog, then imple- 
ment the same background-correction algorithm used on 
the SDSS data to characterize the galaxy population sta- 
tistically associated the MaxBCG-identified cluster cen- 
ters. In order to test the algorithms, we compare the 
properties of galaxies located in three dimensions around 
MaxBCG-identified cluster centers with the statistical 
distribution of galaxies determined by the background- 



4. RESULTS 

Our primary results comprise the luminosity and color 
distributions for cluster galaxies as a function of clus- 
ter richness, cluster-centric distance, and redshift. We 
examine BCG and satellite galaxies separately, moti- 
vated by the observational and theoretical reasons pre- 
viously discussed. As an example of BCG-satellite dif- 
ferences. Figure [5] shows the luminosity function for all 
galaxies within r2oo of cluster centers for four different 
bins of cluster richness. The luminosity functions of the 
satellites (filled circles) and BCGs (open diamonds) are 
shown separately; the total LF is given by the black line 
(with error bars omitted for clarity) . Examination of the 

/d.o.f. reveals that the satellite LF is statistically well 
described by a Schechter function, but this behavior is 
not the case when the BCG is included. The BCG-only 
LF is well-fit by a single Gaussian. In lower richness clus- 
ters, the BCGs dominate the total light but are system- 
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Fig. 7. — The luminosity function of all, red, and blue satellite galaxies in clusters (within r20o), as a function of cluster richness. Each 
panel represents a bin in cluster richness as given in the legend. Black filled circles show the LF for all satellites, with the best-fit Schechter 
function to those data shown as the solid black line. The red and blue galaxy samples arc shown by red diamonds and blue squares 
respectively. While the majority of satellites are red, there is a population of blue galaxies in clusters of all richnesses. Also shown is the 
Schechter function with fixed faint-end slope (o(,iue = 1-0) fit to the blue galaxy distribution (solid blue line). For reference, the best-fitting 
Schechter function for the mean LF of all satellite galaxies in clusters where the catalog is known to be complete and pure {N200 > 10) 
is reproduced in each panel for all satellites (dashed black line) and blue satellites (blue dot-dashed line). The Schechter fits to the red 
galaxies are omitted for clarity. 



atically fainter and have larger luminosity variance than 
BCGs in larger systems. These differences illustrate one 
way in which BCGs tend to be distinct from satellites. 

We begin by presenting results describing the satellite 
population, then turn to the BCGs. In § 14. 1[ we exam- 
ine the luminosities and colors of satellites. Specifically, 
we measure the luminosity function of satellies within 
r2oo as a function of cluster richness; this measurement 
is closely related to the conditional luminosity function 
that is parametrized as a function of mass. We also in- 
vestigate the CLF of red and blue satellites separately, 
and show the radial trend of the LF for all, red, and blue 
satellites. In addition, we measure the red fraction of 
satellites and investigate its dependence on several clus- 
ter properties. In N4.2I we explore the BCG population 
and quantify the relationships of BCG luminosity to to- 
tal cluster mass and to satellite total and characteristic 
luminosity. 



4.1. Satellite Cluster Galaxies 

In this section we measure the conditional luminosity 
function and its dependence on galaxy color and cluster- 
centric distance. We also study how the red fraction 
/r depends on redshift, cluster mass, cluster-centric dis- 
tance, and galaxy luminosity. 

To split the galaxies into red and blue subsamples, we 



make a cut in color-luminosity space: 

0-25(.g -r)= Q.17log,„C-''Lj{LQh^)) - 0.6. (11) 

Because of the strong bimodality in the population, our 
results do not depend sensitively on the exact placement 
of the red-blue boundary. Figure [6] shows the bivari- 
ate distribution of color and luminosity for four example 
richness bins and two radial ranges: 3 < ^''200 < 5; 10 < 
N2ao < 15; 30 < N200 < 40; 80 < N200 < 200 and 
0.05 < r/raoo < 0.3 (left set of 4 panels); 3.0 < r/raoo 
< 5.0 (right set of four panels) respectively. The red se- 
quence exists at all richnesses. Even at large radii, there 
is an excess over random of galaxies associated with clus- 
ters, and many of them are blue. The dashed line in the 
figure shows our adopted cut between red and blue galax- 
ies. 

4.1.1. Conditional Luminosity Functions 

We start by investigating the dependence of total, as 
well as red and blue split, luminosity functions on clus- 
ter richness and cluster-centric distance. We use bins 
defined by sharp cuts in richness; due to scatter in the 
mass-observable relationship, each bin therefore contains 
clusters from a non-sharply defined mass range. Note 
that these observationally-defined conditional luminosity 
functions, which depend on both richness and cluster- 
centric distance, are related to but differ somewhat in 
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Fig. 8. — Parameters of the best-fitting Schcchtcr functions 
from Figure [7] for all satellites (left) and for red and blue satel- 
lites (right) as a function of cluster richness. Symbols are the 
same as in Figure [7] Top: Characteristic galaxy luminosity, L» ; 
middle: faint end slope, a, which is fixed to 1.0 for blue satellites; 
bottom: normalization of the LF, 0*. The best-fitting power law 
describing the Schechter parameters as a function of richness for 
all satellites is shown as a solid line where the fit was performed, 
and as a dashed line extending the trend to lower richness. The 
power law parameters may be found in Table [T] 

definition from the traditional definition of the CLF as 
a function of cluster mass wher e the binning is d one in 
specific mass ranges (as in, e.o.. lYang et al.ll2003[ ). 

Figure [7] shows the luminosity function of satellite 
galaxies within r2oo for different bins of cluster richness. 
These LFs are shown for all satellites (filled circles) and 
for the blue and red subsamples separately (open squares 
and open diamonds respectively) . The shape of the over- 
all satellite LF is only weakly dependent on A''2oo for the 
top two rows, representing higher N200 clusters. These 
arc the systems that have been well studied and are 
known to be quite complete and pure. For A^200 < lOi 
sub- and supcr-L* satellites arc both found with lower 
density than in more massive systems, but the incidence 
of L* galaxies is about the same as in the higher-mass 
objects. With the split into red and blue subsamples, 
we see that change in the faint end of the overall LF is 
driven by the changing ratio of red and blue galaxies. 
We investigate the changing fraction of red galaxies as a 
function of galaxy luminosity and cluster mass further in 

§11321 

We fit the above LFs with lSchechterl ()1976l ) functions, 
(t){M)dM = 

O.41n(lO)0aO-°-4(^"-*^*)("+i)e-i"""''""*'dAf(12) 

using the Levenburg-Marquardt minimization proce- 
dure. The solid black lines in Figure [7] show the best- 
fitting Schcchtcr function for the total satellite LF in each 



TABLE 1 

Power Law Fits to Schechter Parameters 



Parameter 


Units 


Vs. 


N ormalization 


Index 




1O"'L0 


N200 


0.8 ± 0.1 


0.15 ± 0.04 




lOi«L0 


Mi4 


1.29 ± 0.02 


0.12 ± 0.03 


a 




-0.28 ± 0.06 


0.25 ± 0.07 






Ml4 


-0.63 ± 0.02 


0.20 ± 0.05 


0. 


Mpc-3 


^''200 


8 ± 1 


-0.20 ± 0.04 




Mpc-3 


Mi4 


4.3 ± 0.1 


-0.16 ± 0.03 



Note. — Mm = M2oo/(1O"M0) 



richness bin. The satellite LF is well fit by a Schechter 
function except for in the very lowest richness bin. For 
reference, the best-fit Schechter function to the mean 
satellite LF for clusters with A^200 > 10 is repeated in 
all panels as a dashed line. We also fit a Schechter 
function to the LF of the blue and red galaxy popu- 
lations. For blue satellites, in the cases where the fit 
is well constrained, we find the faint end slope to be 
Oiblue = 1-0. In some cases the blue data do not permit 
a strong constraint on the faint-end slope. We therefore 
fit the blue LFs with a faint-end slope fixed to be 1.0. 
These Schechter functions for blue satellites is shown in 
each panel as the solid blue line, and for comparison, the 
best fit for the mean blue satellite LF of clusters with 
-^200 > 10 is repeated in all panels as a dot-dashed line. 
The fits for the red satellites are omitted for the sake of 
clarity. 

The parameters of these best-fitting Schechter func- 
tions are shown in Figure [Has a function of -/V2oo- For 
systems with 7V200 > 10, for all satellites is a weak 

function of mass, and is roughly 1.4 x \Q^^h~^LQ. For 
blue galaxies in clusters of any richness, is con- 

sistent with 1.4 X lO^°/i~^-L0; red galaxies have con- 
sistent with this value only for 7V200 > 30 {W^h-^Mo), 
while lower mass systems have red satellites with a fainter 
characteristic luminosity. The faint end slope, a, scales 
with mass for all and red satellites, with a steeper faint 
end slope for more massive clusters. For blue galaxies, 
we have fixed the value of a, but by inspection of Figure 
[7| it is seen that auue changes little over the whole rich- 
ness range of the catalog. The normalization of the LF, 
0*, decreases with mass for all, red, and blue galaxies. 
For each of L*, a and (/)*, we fit the trend of the param- 
eter as a function of richness for 7V200 > 10. The fits are 
shown on the Figure, with a solid line where the fit was 
performed and a dashed line extending the relationship 
to lower richness. The parameters of these fits are listed 
in Table [H 

These CLF results show the weak dependence of the 
LF parameters on cluster mass, and are in reasonable 
agreement with other recent modeling and observational 
results. The relationship between characteristic satellite 
lu minosity and c luster richness was recently discussed 
b v lSkibba et al.l (I2007T) . They explored the predictions 
of iSkibba et al.l (|2006[ ) that, for halos with mass greater 
than lO^^/i~^M0, their Halo Occupation Distribution 
(HOD) models predict that the mean luminosity of non- 
central galaxies should be nearly independent of mass 
and that the shape of the LF is a pproximately indepen- 
dent of mass. ISkibba et al.l lj2007t ) showed that the mean 
non-central galaxy luminosity is indeed only very weakly 
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Fig. 9. — The luminosity function of all (black filled circles), red (red open diamonds), and blue (blue open squares) satellites as a 
function of distance from cluster center (increasing across rows) and of cluster richness (decreasing down columns). The best-fit Schechter 
function to all satellites in the innermost radial bin of each richness range (solid black lines in left-most column) is rescaled by 0* for the 
other radial bins (dashed lines). The excellent match between this scaled LF and the data illustrates that the LF changes normalization, 
but not shape, as a function of radius for these magnitudes. However, the fraction of red galaxies steadily increases with increasing distance 
from the center, most noticeably for sub-L« galaxies. 

dependent on cluster richness (as quantified by the num- 
b er of galaxie s with Mr < —19.9) in the grou p catalogs 
of lYang et all (l2005l) andlBerhnd et all pOOl ). The cat- 
alogs used in iSkibba et al.l (|2007f ) just barely overlap in 
mass with the lower mass end of our ^ 10 sample. 
Our findings for ^ 10 ai'e in agreement with the 

group-scale results, and show that the HOD prediction 
continues to hold for greater mass s ystems a s well. The 
results of CLF modeling discussed in ICooravl (|2005[ ) and 
iCooravl l|2006[ ). based on constraints from the two-point 
correlation function, also qualitative agree with this the 
behavior of as a function of cluster mass. 

The dependence of the fain t end slope on iV2o o has 
also been discussed elsewhere; ICooravl ()2005l 120061 ) also 
found a steeper faint end slope for more massive sys- 
tems, and a changing f aint end slope was incorpo rated 
into the CLF models of Ivan den Bosch et al.l ()2005l ) in a 
similar way to fit the observed trends in 2dFGRS data. 
The trend for the faint end falling off more rapidly in 
less massive system s is qualitatively sim ilar to the trend 
seen m RCS-1 data (jCilbank et al.|[2007[ ). although more 
work is required to make a detailed comparison with this 
higher-redshift, and very differently selected, sample. We 
note that the magnitude limit of our data is a relatively 
bright "-^^Mi - 51ogio/i < —19, so we cannot comment on 



the faint end turn-u p seen in deeper studies (lYagi et al] 
120021: iPopesso et all 120061: iBarkhouse eTall I2007D . Our 
results for all, red, and blue satellites are in good agree- 
ment with these deeper results in the magnitude range 
where we overlap. 

The somewhat counter-intuitive trend of decreasing 0, 
in more massive systems reflects the correlation between 
(M/L)200: the mass-to-light ratio measured at r2oo, and 
cluster mass. The (M/L)2oo--^20o relatio nship is both 
expected theoretically in some models f e.g.. [Tinker et all 
2005; van den Bosch et al. 2 007., and referen ces therein) 
and directly observed (jSheldon et al.|[2"007bD . By defini- 
tion, the mass density within r2oo is the same for clusters 
of all masses. Since (M/L)2oo increases as a function of 
cluster mass, the luminosity density must be decreasing 
with mass. As the typical galaxy luminosity increases 
with mass, the mean number density of galaxies within 
r2oo must decrease with mass. 

In I Hansen et al.l ()2005D we explored the LFs of galaxies 
in clusters as a function of richness, and it was not clear 
that all satellite LFs within r2oo were well described by 
a Schechter function, especially for the low-richness sys- 
tems. However, the cluster finder used in that previous 
work used a fixed 1 Mpc aperture for estimating richness, 
resulting in much greater scatter in the mass-richness re- 



14 



Hansen et al. 



1.6 - 



C 

CO 

0.8 



0.6 

0.4 



All 




70 II Naoo 


II 200 


40 II Njoo 


II 69 




II 3P 


20 II Ngoo 


II 29 


in II N,^„ 


II 19 


6 II Ngoo II 


g 


3 II ti^oo II 


5 




I nil-,.. 

Blue Sample 



0.1 



1.0 



0.1 



1.0 



r/r 



200 



Fig. 10. — Mean "'^^(g — r) color as a function of r/r200 changes due to the changing red and blue fractions. Left: all satellites; right: 
the sample split into red and blue galaxies. The mean color overall changes, but as the mean color of the red and blue samples does not 
change, the overall shift must be due to the changing numbers of red and blue galaxies. The estimated global field value is shown by the 
dashed line. 

ies to the total LF at large radius is driven by changes in 
the faint end slope of the red LF. The shape of the blue 
galaxy luminosity function changes relatively little with 
radius, and at all radii the blue LFs show proportion- 
ately more faint galaxies than do the red or total LFs. 
The relative constancy of the bright end of the LF of 
all satellites as a function of radius is in agreement with 
the recent detailed examina t ion o f 57 low-redshift Abell 
clusters bv lBarkhouse et all ()2007l ). as is the finding that 
the LFs of red galaxies change more dramatically with 
radius than do the blue LFs. 

These trends are discussed further in the context of the 
red fraction measurements below. 

4.1.2. The Red Fraction & Other Color Statistics 

From examination of the LFs, it is clear that the red 
fraction changes with a variety of parameters. In this 
section we examine the dependence of the red fraction on 
cluster richness, cluster redshift, cluster-centric distance, 
and satellite galaxy luminosity. 

To assess whether changes in the galaxy population are 
primarily related to changes in the red fraction, we exam- 
ine the typical color of satellites as a function of distance 
from the BCG. Figure [TOl shows the (number weighted) 



lation than in the present version of the cluster-finding 
algorithm, thereby mixing together clusters of much dif- 
ferent mass into the same bin of richness. In addition, as 
we will show, the LF normalization changes dramatically 
as a function of r/r2Qo, so a fixed aperture will mix popu- 
lations from different fractional r/r2oo sections of clusters 
and thus further confuse the analysis. While any catalog 
will have noise and impurities, the current cluster cata- 
log, which uses and improved richness measure, should 
be less affected by this issue. 

In addition to the LF within r2oo! we also investigate 
the radial dependence of the luminosity function for all 
galaxies and for red and blue galaxies separately in sev- 
eral bins of cluster richness. These results are shown in 
Figure m Cluster richness decreases downward in the ar- 
ray of plots, and the radius range increases to the right. 
There are still galaxies correlated with the clusters cen- 
ters at separations much larger than r2oo- While these 
galaxies are likely not bound to the clusters, they are 
statistically correlated with clusters and the amplitude 
of the luminosity function reflects the amplitude of the 
cluster-galaxy cross-correlation function. 

The normalization of the LF decreases with cluster- 
centric distance, but the shape of the total satellite LF 
remains essentially unchanged for the magnitude range 
considered here. To illustrate this point, we have re- 
drawn the best-fit Schechter function for all satellites in 
the innermost radial bin (solid lines in left hand column) 
in the panels for other radial bins, but rescaled by (p^ 
(dashed lines in other columns). This rescaled LF is an 
excellent match to the measurements. For reference, we 
find that in the 0.25 < r/r2oo < 0.75, 0.75 < r/r2oo < 
1.5 and 1.5 < r/r2oo < 5 radial bins respectively, is 
typically 20%, 3%, and 0.3% of in the innermost bin. 

Although the relative mix of luminosities is unchanged 
with radius, there is a strong change in the relative mix of 
colors. Red galaxies are much more dominant in the in- 
ner regions, and the decreasing contribution of red galax- 



mean {g — r) color for satellite galaxies with ' °A-I, 
51ogio/i < —19, in bins of cluster-centric distance out to 
5xr2oo, and for several bins in cluster richness. The left 
panel shows the trend for all galaxies, while the right- 
hand panel shows the mean color for the red and blue 
samples of galaxies separately. The dashed line shows 
the estimated mean color of the universal average popu- 
lation. The overall mean satellite color is redder in the 
inner regions and trends bluer as a function of radius 
until r ~ 2xr2oo, beyond which the mean color is con- 
sistent with the universal average value. The mean color 
is insensitive to cluster richness at any radius, with the 
exception of the very lowest richness bin in which the 
galaxies are typically bluer than in richer clusters. How- 
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ever, the right panel shows that the red and blue samples 
remain the same color regardless of cluster-centric dis- 
tance or cluster richness, which implies that the overall 
shift in color is simply a change in the relative number of 
red and blue galaxies. This observation suggests either 
that there is some process that distributes red and blue 
galaxies differently with radius, or that the transition 
between the two populations happens fairly rapidly. Re- 
gardless of the mechanism, the color shift appears to be 
characterized well by the changing number of red galax- 
ies, and we subsequently consider the red fraction, /r , as 
the primary statistic for quantifying the change in galaxy 
population. 

We first measure the satellite red fraction within r2oo 
as a function of cluster richness and redshift. Figure [Til 
shows the fraction of red satellites brighter than °-^^A/.j 
- 51ogioft. = — 19 and within r2oo as a function of clus- 
ter richness. The clusters are split into two bins of red- 
shift: 0.1 < z < 0.25 (filled circles, median z = 0.2) 
and 0.25 < z < 0.3 (open diamonds, median z = 0.28). 
The red fraction increases significantly as a function of 
cluster mass in both redshift slices, but for lower redshift 
systems /r is systematically higher at all richnesses by 
~5%. There is a hint that /r, flattens above W^hr^MQ. 
To parametrize the trend of red satellite fraction as a 
function of cluster richness and redshift, we adopt the 
functional form 

!r = g{z) erf [logio(A^20o) - \o^Mm + 0.69. (13) 

We find that with 

g(z) = 0.206 - 0.371Z (14) 

/i(z) = -3.6 + 25.8Z (15) 

this form provides a reasonable description of the data; 
this relationship is shown as the solid lines in Figure [TTl 
For A^200 < 10, the catalog may be less than 90% com- 
plete and pure, and it is possible that as-yet-unquantified 
selection effects may influence the trend in this regime. 

Here, we have made measurements within the aperture 
r2oo measured for the full cluster sample with the lensing 
data, relative to the critial density at z = 0.25. However, 
as the critical density scales with redshift, one can ques- 
tion whether we have biased our result by not using the 
r2oo taken relative to the critical density appropriate for 
the mean redshifts of the two subsamples. While the 
lensing analysis has not been done separately for the two 
redshift ranges, we can investigate the dependence of r2oo 
on redshift using the galaxy density estimate. We find 
that using the z-appropriate threshhold for the high and 
low redshift samples results in r2oo values that are 8% 
and < 1% different, respectively, from the value found 
for the whole sample. Using an adjusted aperture size, 
we measure red fractions that are 0.5% different for the 
higher redshift sample. The change in the low redshift 
sample is even smaller. This systematic offset is clearly 
not sufficient to explain the observed 5% difference in red 
fraction between the the two redshift bins. 

As seen from the behavior of the red and blue LFs, the 
red fraction has additional dependence on cluster-centric 
distance and satellite galaxy luminosity. In Figure [T^] 
we examine both the radial and luminosity dependence 
of the red fraction. The left panel of Figure [T^] shows 
/r(''A20o) for satelUies with "-^^M; - 51ogio/i < -19 for 
clusters in several bins in richness. Within ^ 2xr2oo, 
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Fig. 11. — Red fraction within r2oo (°'^^Mi - Slogio/i < —19) as 
a function of richness for two redshift bins: 0.1 < z < 0.25 (black 
filled circles, median z = 0.2) and 0.25 ^ z < 0.3 (brown open 
diamonds, median z = 0.28). The solid lines arc the richness and 
redshift dependent model of Equation 1131 

the red fraction decreases with radius. Beyond 2xr2oo, 
/r flattens asymptotically, approaching the cosmological 
average (dashed line). The red fraction is essentially in- 
dependent of cluster richness for A^200 > 10 systems, but 
is correlated with richness for 7V200 < 10. 

The right panel of the Figure shows /r for satellites 
within r2oo as a function of absolute magnitude for multi- 
ple richness bins. For satellites brighter than L*, the red 
fraction is relatively independent of luminosity except for 
the lowest richness systems. For sub-L* galaxies, the red 
fraction steadily decreases toward fainter magnitudes, 
and the trend is stronger for lower iV2oo systems. Ex- 
cept for the brightest galaxies in the very lowest-richness 
bin, /r is always greater than the universal average value 
(dashed line). 

Broadly speaking, the fraction of red satellite galaxies 
increases with cluster mass, with galaxy luminosity, with 
time, and for decreasing cluster-centric distance. The 
dependence on cluster richness is coupled more strongly 
to satellite luminosity than to cluster-centric distance. 
The previous LF results showed that the faint ends of 
the red and blue satellites' LFs change as a function of 
cluster-centric distance, although the bright ends remain 
essentially unchanged. In light of this observation, we 
conclude that it is the steadily changing mix of sub-L* 
galaxies as a function of distance from the cluster cen- 
ter that causes the trend of decreasing /r with cluster- 
centric distance. 

These results are in general agreement with pre- 
vious work, despite differences in sample definitions 
for both cluster selection and sample sp lit ting. Sev- 
eral studies (iMargoniner et al ] l2001t iGoto et all 
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|2007| ) have reported that the galaxy type fraction is a 
function of cluster mass. How ever, others (iBalogh et al.1 
I2004bt iDe Propris et al.ll2004l : iTanaka et al.ll2004h found 



Fig. 12. — Dependence of red fraction on cluster-centric distance and galaxy luminosity. Left: Red fraction (for galaxies with " ■^^M; - 
51ogio/i < —19) as a function of r/r2oo for different richness bins. The estimated global field value is shown by the dashed line. Right: Red 
fraction (within r2oo) ^ function of absolute magnitude for different richness bins. The estimated global field value for each luminosity 
is shown by the dashed line. 

forward because of the very different sample selection 
criteria used, but in future work we will undertake this 
comparison. 

4.2. Luminosity of Brightest Cluster Galaxies 

We now examine the Brightest Cluster Galaxies, fo- 
cusing on their luminosities and how those luminosities 
compare to the total cluster luminosity and the charac- 
teristic satellite luminosity. 

The top panel of Figure [13] shows the median luminos- 
ity of BCGs, °-^'^LbcGj as a function of cluster richness. 
Systems with a greater number of satellites tend to host 
brighter BCGs. The error bars on the data points indi- 
cate the statistical uncertainty on the median BCG lumi- 
nosity in each bin, determined from jackknife resampling; 
the dashed lines show the region within which 68% of the 
BCG luminosities lie. For A^200 > 10, we parametrize the 
^'^^LscG-richness trend with a power law, and find that 
the z-band light from BCGs scales with cluster richness 
as 

°-^'Lbcg = (2.16 ±0.08) X 10i"/i-2LQiV2Vo'" °'. (16) 

This fit is shown as a solid line Figure 1131 Using our 
adopted mass-observable scaling, this relationship is 

0.30±0.01 



that the type fraction does not depend on cluster 
velocity dispersion, although iWeinmann et al.l (|2006aD 
argues that this result may be due to the large scatter 
in relationship between velocity dispersion and mass. 
Furthermore, the flattening and "saturating" of the 
red fraction for the most massive systems appears to 
be a natural prod uct of modeling (iBerlind et al.l l2003l: 
Zheng ct al] |2005|). Some a uthors (iDc Propri s et al.l 
2004 : ".Weinmann et all l2006al : iBlanton fc Berhndl l2007f l 
have also found qualitatively similar behavior for the 
type fraction as a function of radius; however, this 
trend is very weak in both thc lWilman et all ()2005f ) and 
iGerke et al.l (|2007f ) data, but there are large uncertainties 
in both of those data sets that could mask a correlation. 
The luminosity dependence of the gal axy type fraction 
is al so similar to p r evious results (D e Propris et all 
2004t Iwilman et al'l l2005l: IWc inmann et all l2006al: 



Martinez et al.ll2006l : ICerke et al.li2007t ) 



The redshift dependence of the gala xy type fraction has 
been d iscussed by many authors since iButcher fc Oemlerl 
([1981 . OveraU, we find that the fraction of red galaxies 
increases by ~5% during the 0.8 Gyr between our higher 
(median z = 0.28) and lower (median z = 0.2) redshift 
samples over the full mass range. This result is in quite 
good agreement with the results of the eponymous work, 
where an increase in the blue fraction by 0.25 in clusters 
to z = 0.5 was observed ([But cher & Ocmlcr 1984). We 
also find agreement with [Goto et al.l (|2003|), who used a 
smaller catalog of clusters from the SDSS EDR. With 
the larger MaxBCG sample, however, we can see that 
this shift in red fraction happens in a similar manner for 
clusters of all richnesses. Direct comparison with higher 
redshift samples, suc h as the blue fracti on of the DEEP2 
groups measured bv lGerke et al.l (|2007f ). is not straight- 



° ''ii3CG = (6 ± 2) X m'"h'^LQ 



200 



Ml 



(17) 

where we ha ve defined M14 = \Q^ ' ^ h~^M(r,. A scaling 
suggested bv lVale fc Ostrike"i1 (|200(3D for the LBCG-mass 
relation is 



(18) 



'[l-f (M/Mc)'''=]i/fe' 

with A4 = 3.7 X IO^/i-^Mq, a = 29.78, b = 29.5, and 
k = 0.0255. This functional form is derived by fitting 
the relationship obtained by matching the galaxy lumi- 
nosity function to the subhalo mass function, and is not 



The Galaxy Content of SDSS Clusters and Groups 



17 



5.0 10.0 



10 



T, 0.24 

S 0.22 

J° 0.20 

60 0.18 

o 0.16 

^ 0.14 

1= 0.12 




1 



10 



100 



Fig. 13.— Top: Median BCG O-^^i-band luminosity, '' ^SLscg, 
as a function of richness. Error bars show the statistical uncer- 
tainty on the luminosity; the dashed lines show the 68% scatter in 
^'^^ LgcQ in each richness bin. The solid line shows the best-fit 
power law, '^'^^Lbcg ~ Mgn^, fit for Af200 > tO- T he dot-dashed 
line shows the best-fitting IVale &: Ostrikeri ~ l|2006l 'l model. Bot- 
tom: The width of the 68% scatter of ^'^^ Lbcg a function of 
richness. For massive systems this width is ~ 0.17. 

currently physically motivated. Equation [18] also pro- 
vides an acceptable fit to our data if we adopt Lq = 
4 X lO^/i~^L0, the mean galaxy luminosity in halos of 
3.46 X 10"/i-iMq in the " ^^i-band. This relationship is 
shown on the Figure with the dot-dashed line. 

The trend of increasing central galaxy luminosity with 
cluster richness has be en noted in many previous obser- 
vational studies (e.g.. ISandage fc HardvlTl973l: ISandagd 
[TOTfil: iHoessel et alTllQSnl : iSchneider et all I1983D . New 
large cluster samples with robust mass estimators have 
explored the scaling of BCG luminosity with cluster 
mass. Using a sa mple of 93 c l usters with both X-ray 
and if-band data, iLin fc Mohd (|2004f ) found that BCG 
light scales with mass as LBCG,K—band ~ 

AfOo§^ for clus- 
ters with M200 > 3 X 1 0^^/i.~^M(7), with significant scat- 
ter. lYang et all (l2005h . using groups found in the 2dF- 
GRS, found that in &j-band, (L ^ p„.) ~ M°-^^ for halos of 
M > W^h-^MQ. I Zheng et all |2007D . in their investi- 



gation of the luminosity-dependent projected two-point 
correlation function of DEEP2 and SDSS, also found 
that there is a correlation betw een halo mass and cen - 
tral galaxy luminosity. Using the I Vale fc Ostrikeil ()2006f ) 
model, they see that Lq = 2.8 x lO^/i~^L0 provides a 
reaonable fit for the SDSS r-band (with halo masses up 
to 3 X IQI^/i-IMq) and Lq = 4.3 x lO^/i-^ig is suitable 
for the DEEP2 B-band data (with halo masses up to 
4 X \Q^^h~ ^M(i,). Using a differ ent cluster catalog from 
the SDSS, iPopesso et all (|200l found Lbcg M^i^f 
for these 217 systems. Our results are in agreement with 
these findings within the uncertainties, but the size of 
the MaxBCG catalog, its well-understood selection func- 
tion and its accurate mass estimator allow us to probe 
the BCG luminosity distribution in further detail. 



At all richnesses, the distribution of ^'^^Lbcg is well- 
described by a Gaussian. The mean value is dependent 
on cluster richness as previously discussed. The width 
of the distribution, <JiogL, is also a function of cluster 
richness. The bottom panel of Figure [13] shows aiogL 
as a function of richness, with error bars from jackknife 
resampling. There is an overall negative correlation be- 
tween aiogL and iV2oo, although for ~35, aiogh is 
roughly consistent with a constant value of ^ 0.17. This 
val ue is somew h at hig her than the aiogL ~ 0.12 found 
by IZheng et al.l (|2007 ). but is consistent within the un- 
certainties. Our measurements are made as a function 
of cluster richness, and scatter in the mass-observable 
relation causes clusters over some range of masses to be 
assigned to each richness. Since ^'^^Lbcg depends on 
cluster mass, this scatter may result in larger observed 
o'logL values than in the intrinsic distribution of BCGs as 
a function of cluster mass. However, mass mixing acts 
only to increase the observed value, so our measurements 
represent at least a secure upper limit on the scatter in 
the '^'^''^ Lbcg~M2oo relationship. Future work with simu- 
lations is required to fully disentagle the intrinsic scatter 
from that introduced by the mass proxy. 

The MaxBCG cluster finder includes priors on BCG 
color and luminosity (see § 12. 2p . and so we consider 
whether the resulting BCG luminosity distribution is be- 
ing artificially constrained by these priors. In most cases, 
where BCG identification is unambiguous, only the nar- 
row color priors inform BCG selection, so we do not ex- 
pect a significant effect from the magnitude prior. De- 
tailed examination of the effect of these priors on the 
selection function of MaxBCG is underway, but prelim- 
inary results indicate that the incidence of rich clusters 
with ambiguous BCGs is < 25%, and does not have a 
significant effect on the scatter of rich systems. Here, 
we see that the width of the distribution of identified 
BCG i-band absolute magnitudes is ~ O.Smag wider than 
the (mostly uninformativc) magnitude prior, and that all 
BCGs are easily brighter than the nominal 0.4iH. limit. 
Thus, we take the recovered distributions as representa- 
tive of the cluster population, but reserve a more detailed 
investigation for future work. 

In addition to examining the correlation between BCG 
luminosity and cluster mass, we also measure the trend 
with mass of the ratios of BCG luminosity to the to- 
tal cluster luminosity, ^2001 and BCG luminosity to the 
characteristic luminosity, ''■^^i*(sat), of the satellites. 
The top panel of Figure [14] shows ^'^^Lbcg/ as a 
function of cluster richness. As expected from the exam- 
ple LFs of Figure [5] °'^^Lbcg/L2oo decreases with clus- 
ter richness. For the most massive clusters (IQI^/i-IMq), 
the BCGs supply only ~ 5% of the cluster luminosity 
budget, but for intermediate IO^'^H^^Mq systems the 
BCG makes up 20% of the light. For the lowest rich- 
nesses the luminosity is completely dominated by the 
BCG. 

Fitting a simple power law, we find that for clusters 
with iV200 > 10, the BCG Ught fraction scales with clus- 
ter richness as 
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Fig. 14. — BCG luminosity compared to the luminosity of the rest 
of the cluster galaxies within r2oo ("'^^M^ - 51ogio/i < —19). Top: 
Mean BCG luminosity fraction C^'^^ Lbcg/ L200) as a function of 
richness. Bottom: The ratio of BCG luminosity to L* of the 
satellite galaxies, ^'■■^^ LgQQ/^''^'' Lt{sat), as a function of richness. 
In each case the best-fitting power law is shown with a solid line 
where the fit was performed and a dashed line to extend the relation 
to lower richness. 



and with cluster mass as 
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This relationship is shown in the Figure with a solid line 
over the A'^200 range where the fit was performed, and as 
a dashed line to extend the relation to lower richness. 

The negative correlation between the BCG luminosity 
fraction and cluster mass is in good agreeme nt with other 
result s. The K-hand data presented by iLin fc Mohrl 
()2004f l for their X-ray selected sample expresses quali- 
tatively the same tr end. For a different catalog of X- 
ray selected clusters, iPopesso et al.l (|2007[ ) measured the 
observed scaling between total optical luminosity and 



cluster mass and between BCG luminosity and cluster 
mass; the combination of these relationships yields a best 
fit scaling of ^'^^ Lbcg/ L20Q ~ -^^200^^' consistent with 
our results. We are also in reasonable agreement with 
the low redshift { z = .1), r-band model predictions of 
iTasitsiomi et al.l (j2007[ ). who expect ^ "^^ LbcgI ~ 
massive systems, and who also predict that 
the normalization of this relationship may be sensitive to 

(Is- 

We also compare "'^^iscG to the characteristic satel- 
lite luminosity °-^^L*(sat). The value of °-^^L*(sat) is de- 
rived from fits to Schechter functions for satellites within 
r2oo and brighter than °-^^Afi - Slogioft. = —19 as in § 
14.1.11 The bottom panel of Figure [14] shows that the ra- 
tio ^-"^^ Lbcg /^'"^^ L^,{sat) depends strongly on A^200- For 
clusters with A^200 > 10, ^ '^^ LBGG/^'^^L^{sat) increases 
with N20Q, but this trend reverses for very low richness 
systems. This behavior is not unexpected given the re- 
sults of Figures [8l and [T3l '^'^^Lbcg nionotonically in- 
creases with iV2oo, while '^■'^^^ L^,{sat) is essentially flat for 
-^200 > 10 but correlated with richness for lower iV2oo 
systems. We caution that the low richness (iV2oo < 10) 
systems should not be taken as representative of all low- 
mass systems, as the selection function imposed by the 
cluster finder (demanding only a few red galaxies in close 
proximity) may be important in this regime. 

For 7V200 > 10, the ratio °'^^LBCG/°'^^L^{sat) as a 
function of cluster richness can be described by a power 
law: 
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and likewise as a function of mass: 
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This relationship is shown in the Figure with a solid line 
over the range in N2oa where the fit was performed and 
a dashed line to extrapolate the trend to lower richness. 

5. SUMMARY & DICUSSION 

In this study, we have examined the properties of 
cluster-associated galaxies, separating BCGs from satel- 
lites and focusing on trends in galaxy color and lumi- 
nosity as a function of cluster richness, distance from 
cluster center, and redshift. We use clusters in the SDSS 
MaxBCG sample, the largest set of clusters identified to 
date. Employing photometric data alone, we apply cross- 
correlation background-correction techniques to charac- 
terize the cluster-associated galaxy population, within 
5xr2oo and brighter than °-^^Mi - 51ogio/i < —19, around 
165,597 systems spanning more than two decades in mass 
in the redshift range 0.1 < z < 0.3. 

Our principle results are as follows. 

1 . The luminosity function of satellites within r2oo as 
a function of cluster mass for systems with mass 
greater than 3 x 10^^h~^ Mq shows remarkable uni- 
formity for °-^^Mi - 51ogioft. < —19. The character- 
istic satellite luminosity L» is only weakly depen- 
dent on cluster richness. 

2. The shape of the luminosity function of satellites 
brighter than "-^^Mi - 51ogio/i < —19 does not 
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change with cluster-centric radius. However, the 
color-separated luminosity functions of satellites as 
a function of r/r2oo and of iV2oo show that the mix 
of sub-L* red and blue galaxies changes dramati- 
cally as a function of radius. In contrast, the rela- 
tive number of red and blue galaxies at the bright 
end is roughly constant with radius. 

3. The average color of satellite galaxies is redder near 
cluster centers, but this trend is largely a reflection 
of the changing ratio of red to blue galaxies men- 
tioned above. This effect is a very weak function 
of cluster mass over the range investigated here. 

4. The fraction of red galaxies increases with cluster 
mass over the full range explored, although only 
weakly for M200 ^ IO^^'/i'^Mq. This fraction de- 
creases with cluster-centric distance until 2xr2oo! 
decreases with luminosity for L < L*, and is con- 
stant for r > 2x^200 and L > L*. 

5. The fraction of cluster galaxies within r2oo that are 
red increases by ~5% during the 0.8 Gyr between 
redshift z = 0.28 to z = 0.2; this change is roughly 
independent of mass over the range investigated. 

6. The luminosity of BCGs and the ratios of BCG lu- 
minosity to total cluster luminosity and to char- 
acteristic satellite luminosity are all correlated 
strongly with cluster mass, and we have quanti- 
fied each of these scalings over the mass range 
3xl0i3/i-iMo to9xlO"/i-iM0. The BCG lumi- 
nosity has a Gaussian distribution at fixed cluster 
richness, with dispersion oiogh ~ 0.17 for clusters 
with M200 > \^^^h-^MQ. 

While these results are in general agreement with pre- 
vious observational work, due to the volume probed we 
are able to investigate a wider mass range, extending the 
statistics to higher mass than previous samples. We are 
also able to split the sample into finer bins for several 
variables than has previously been possible. 

The MaxBCG selection function, which is well- 
understood for most of the richness range used here, is 
less well quantified for clusters with A^200 < 10- This 
low richness set of clusters may be less complete and 
pure, which could result in the significant difference in 
LF shape as compared to higher iV2oo systems, or con- 
tribute to the drop-off in /r for low richness systems. 
However, not all cluster galaxy properties change signifi- 
cantly at = 10 (e.(?., ^ "^^LbcgI L2m) and there is no 
break at a p articular iV2oo in either the mass-observable 
relat ionship (iJohnston et al.ll2007H l or mass-to-light ra- 
tios (ISheldon et al.ll2007H ) as measured by lensing. In- 
terestingly, it is the quantities that most closely trace 
the total mass of the systems {i.e., ^'^^Lbcg, i20oand 
the lensing signal) that are smoothly scaling over the full 
richness range, while quantities that are related to the 
mix of galaxies within clusters (i.e., the LF and /r) are 
the ones that change more dramatically. We hypothesize 
that, at these low richnesses, MaxBCG is finding legit- 
imate low-mass systems, but that the selection priors 
demanding the close proximity of only a few red galaxies 
result in finding only systems with the observed mix of 
galaxies and not necessarily all systems of this low mass. 



Further investigation using lower mass threshold mock 
catalogs is needed to understand in detail the selection 
of low-richness systems. 

Our results can shed light both on the processes that 
build up the galaxy population in clusters and distinguish 
central galaxies from satellites, as well as the processes 
that are responsible for the galaxy transformation from 
blue to red. With respect to the former, the results pre- 
sented here fit well within the basic picture of galaxy for- 
mation in CDM: that galaxy properties are likely linked 
to the formation history of a cluster's dark matter halo 
and its substructures. Although detailed comparisons 
are beyond the scope of this work, our results qual- 
itatively match both HOD constraints from clustering 
statistics as well as models based on matching the abun- 
dance of halos and subhalos to galaxies. The HOD frame- 
work provides a way to examine the galaxy population in 
both the observed Universe and in models of galaxy for- 
mation and evolution. Without needing to identify spe- 
cific groups or clusters in the data, halo model interpre- 
tations of the statistics of luminosity-dependent galaxy 
clustering result in specific predictions for trends of both 
central and satellite galaxies as a function of halo mass 
that are in reasonable agreement with the findings pre- 
sented here, especially for the relati onship between cen- 
tral and satellite l uminosities {e.g., iBerlind et al .1120031 : 
ISkibba et al.ll2007l ). Models which link the properties of 
galaxies directly to their halos and subhalos also agree 
broadly with several of the results prese nted here {e.g., 
IConrov et ai1l2006l : IVale fc Ostrikedl2007f ) . Detailed pre- 
dictions for several cluster statistics from such a model 
will be presented in lTasitsiomi et all ()2007[ ). Our results 
on scatter in the BCG luminosity at fixed cluster richness 
very likely provide an upper limit on scatter in central 
galaxy luminosity at fixed mass, as the scatter between 
halo mass and cluster richness should act to increase this 
scatter. The fact that this scatter is already fairly small 
provides further support for the tight coupling between 
halo mass and galaxy luminosity that is the basis of these 
models. 

In addition to processes that cause physical changes 
to satellites resulting from interactions with the clus- 
ter gas, cluster potential, or other satellites, presum- 
ably some of the satellites are lost due to being accreted 
onto the BCG. In detail, this process likely results in the 
stellar component of the disrupted ga laxy joining both 
the BCG and ICL (jConrov et al.ll2007f l. but nonetheless 
should result in a BCG population closely linked to both 
halo mass and satellite population, as is observed. In- 
deed, °-^^Lbcg/L2oo and °-^^LBCG/°-^^Msat) must be 
intimately related to the processes responsible for BCG 
growth. That BCGs get brighter as a function of cluster 
mass faster than do typical satellies may be further evi- 
dence that BCGs are different (in their merger history) 
than typical satellites. 

Understanding the timescales and mechanisms for 
galaxies to transform from star-forming galaxies onto the 
red sequence is one of the primary current challenges for 
galaxy formation theories. There are several processes 
that can operate within clusters to shape the population 
of the cluster galaxies, such as ram pressure stripping, 
harassment and strangulation, that directly influence the 
galaxies' gas content and thus t heir subsequent star for- 
mation (for a recent review, see lDe Luciall2006l ). Clearly 
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many of the processes responsible for this transforma- 
tion are related either to the mass of the host halo or 
how long the galaxies have been satellites. The rela- 
tive strengths of various effects are still rather uncertain, 
however. What fraction of galaxies become red while 
they are central galaxies? Do the processes happen only 
for satellites in a certain mass range, or do they happen 
equally for all satellite galaxies? On what timescales do 
these processes operate? Our measurements of how the 
red fraction scales with cluster mass, radius and redshift 
will be instrumental in answering these questions. 

We address a few of these issues here. Ram pres- 
sure stripping predicts that /r, will be inversely corre- 
lated with cluster-centric distance, and larger for both 
brighter galaxies and more massive halos. However, 
our results indicate that is essentially independent 
of galaxy luminosity at fixed iV2oo for galaxies brighter 
than and that the radial trend in /f> is not any more 
pronounced in high iV2oo systems. These results indi- 
cate that ram pressure stripping cannot be the dominant 
process at work to transform the galaxy population. The 
harassment scenario predicts that /r, will be anticorre- 
lated with cluster mass, as is observed; however, if this 
mechanism were dominant, then at fixed cluster mass /r 
wo uld be expected t o be larger for less luminous galax- 
ies ([Weinmann et al . 2006a), contradicting the observed 
trends. Strangulation, where star formation in infalling 
galaxies is halted because no further gas accretion is al- 
lowed, makes several predictions that are in good agree- 
ment with our observations. For example, the model pre- 
sented by Ipiaferio et al. (2 00ll ) predicts: that the mean 
satellite color gets bluer as a function of cluster-centric 
distance until reaching a plateau at the field value around 
2-3r2oo; that the mean satellite color depends on halo 
mass only for M < 5 x IO^^H'^Mq; and that the incidence 
of blue galaxies in clusters increases at higher redshift. 
However, this model also predicts that both bulge- and 
disk-dominated satellies will get redder toward the cen- 
tral regions of clusters, a trend that we do not see in our 
red/blue sample split (or course, our color-separated sub- 
samples are not directly comparable to their morphology- 
separated samples). Note that the predictions of such a 
model will in detail depend on its implementation. Our 
observation that the satellite red fraction changes in the 
same way with redshift regardless of cluster richness is 
significant evidence that the timescale of the physics re- 
sponsible for quenching is the same in systems over the 
full mass range examined. The difference in cosmic time 
between the median redshifts of our two samples is ap- 
proximately the dynamical timescale, over which we ob- 
serve a ~5% change in /r, and this observation should 
allow for useful constraints in studies on the details of 
quenching mechanisms. 

A consensus is emerging from a variety of modeling 
efforts that star formation proceeds most efficiently in 
a mass range around i*, and is less efficient in more 
massive halos and for satellite galaxies. Simple mod- 
els based on this type of assumption can produce many 
of the rough trends that we have seen here; for exam- 
ple, that the mean color of galaxies will not change sig- 
nificantly as a function of halo mass for the range o f 
masses we investigate here {e.g., iDiaferio et al.l [200ll ) . 
and that the galaxy type fraction should b e a weak func- 
tion of halo mass for M >1013/i-1Mq ()Berlind et all 



[2003t iZheng et all [20051 : iCoorav fc MilosavlievidlMsh . 
Furthermore, the general concept of galaxies falling in, 
quenching, and fading matches well with the observed 
radial trends of LF and red fraction. However, most 
of the detailed semi-analytic modeling efforts have had 
some trouble matching in detail observations of the color 
distribution of gal axies and how it changes with envi- 
ronment (see e.g., iCoil et al.l[2007l) . Unresolved issues 
include the rates at which star formation is triggered or 
shut off, and accordingly red galaxies tend to be overpro- 
duced. To understand in detail the physical mechanisms 
responsible for the quenching of star formation as clus- 
ters are assembled, further work is needed to accurately 
reproduce the observed trends in the cluster galaxy pop- 
ulation. The present findings set a local-universe target 
for modeling results, and provide some guidance for the 
relative importance of some of the germane effects. 

A substantial effort over the next decade will be de- 
voted to large-scale, multi-band photometric surveys, in- 
cluding DBS, Pan-STARRS, and LSST. Although the 
primary science driver of many of these projects is to in- 
vestigate the nature of dark energy, the resulting data 
are likely to also provide strong constraints on the pro- 
cesses of galaxy evolution. The results presented here 
provide a low-redshift baseline against which current and 
future high-redshift samples may be compared. From a 
technical standpoint, these data are informative to next- 
generation cluster-finding techniques, and are useful in- 
put for creating the mock catalogs necessary for inter- 
preting cluster surveys. Furthermore, the techniques 
presented here, which use photometric data alone, are 
directly applicable to these upcoming imaging surveys, 
and will thus enable detailed studies of the galaxy popu- 
lation at significantly higher redshifts without extensive 
spectroscopy. 
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